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The detailed dynamics around intermittency bursts is investigated in turbulent shell models. We 
observe that the amplitude of the high wave number velocity modes vanishes before each burst, 
meaning that the fixed point in zero and not the Kolmogorov fixed point determines the intermit- 
tency. The phases of the field organize during the burst, and after a burst the field oscillates back 
to the laminar level. We explain this behavior from the variations in the values of the dissipation 
and the advection around the zero fixed point. 



One of the most fundamental problems in turbulence research is the understanding of intermittency effects [Q. 
In fully developed turbulent flows laminar quiescent periods are interrupted by strong intermittency bursts of high 
energy dissipation. It is well-known both from a number of experiments and numerical simulations that intermittency 
effects cause corrections to the classical Kolmogorov theory j|] when the structure functions of the velocity field are 
statistically averaged over space and/or time [[!]. On the other hand very little is known about the particular structure 
of intermittency, as for instance the shape of and the behavior around a specific intermittent burst. It is the purpose 
of this letter to present a detailed investigation of the behavior of the velocity field before and after a burst takes 
place. We observe a consistent picture in which the velocity gradients, over a small scale, always becomes organized 
and vanish just before the energy burst sets in by an "explosion" in the field. This is like a "calm before the storm" 
and one can draw an analogy to depinning charge-density-waves which phase-organizes just before they slip Q| or 
self-organized-critical systems in general where bursts or "avalanches" of all sizes can be triggered even by the slightest 
perturbation [Q. We draw our conclusions from investigations of shell models for turbulence, which are completely 
deterministic systems where the intriguing structure of intermittency is created by the internal chaotic dynamics. 
Our results indicate that the trivial fixed point in zero and not the Kolmogorov fixed point is responsible for 
the intermittency. Shell models are formed by various truncation techniques of the Navier-Stokes equations and have 
become paradigm models for the study of turbulence at very high Reynolds numbers Q . The mostly studied shell 
model is the "GOY" model of Glcdzer-Ohkitani-Yamada This model yields corrections to the Kolmogorov 

theory || in good agreement with experiments |l5|-|l7|]. 

For the "GOY" shell model, wave-number space is divided into N separated shells each characterized by a wave- 
number k n = A™ fco (A = 2), with n = 1, • • - N, The corresponding amplitude of the velocity field at shell n is a 
complex variable u n . By assuming interactions among nearest and next nearest neighbour shells and phase space 
volume conservation one arrives at the following the evolution equations ^] 

^ t +vk n) U « =*fcn(°n<+l<+2 + y<-l<+l + y K-lK-2) + /<*n,4, (1) 

with boundary conditions b\ — = c\ — C2 = (in-\ = cln = 0. / is an external, constant forcing, here on the forth 
mode. 

The coefficients of the non-linear terms must follow the relation a n + b n+ \ + c n +2 = in order to satisfy the 
conservation of energy, E — J2 n \ u n\ 2 , when / = v = 0. The constraints still leave a free parameter 6 so that one can 
set a n = 1, b n = —S, c n — —(1 — 8) [ p^[ . If helicity conservation is also demanded, one obtains the canonical value 
(5 = 1/2 |l2). The set (Q) of N coupled ordinary differential equations can be numerically integrated by standard 
techniques. In the simulations, we use the following values: S = 1/2, N = 19, v = 10~ 6 , k = 2~ 5 , / = (1 + i) * 0.005. 

Taking a closer look at the dynamics of the GOY model in terms of the complex field u n (t) = r n (t) e l9n w , the 
intermittent bursts consist of a collection of different organizations of the amplitudes r n which travel with exponentially 
increasing speed from the lower up towards the higher shells where they are damped away by viscosity Jl8| . Every 
burst in the model follow a common pattern, where the most prominent characteristic is that the amplitudes of the 
higher modes vanish just before a burst, as shown in Fig.l. During the attraction towards u n = the phases 9 n 
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organize in period three in the shell index n. Just at the point of minimum amplitude, the phases change so that 
a new organization of period three occurs during the rapid repulsion from zero. After the burst the modes oscillate 
back to the laminar level with increasing oscillation periods (Fig.l). 

In order to explain these findings we will use that the period three organization of the phases is present not only 
at bursts, but also during most of the evolution. To argue for this, Fig. 2 shows a long time average of the phase 
differences P n (t) = \9 n — 6* T1 . — 3 1 . We observe, that the assumption holds very well for the highest shells, i.e. that 
6 n -x = 9 n+ 2, #71-2 = #n+i [jn§- With this assumption, the GOY model can be separated into two coupled ODE's 
controlling the evolution of the moduli r n (t) and the phases 9 n (t): 

f n + ir n 9 n = -uk^r n + ik n e~ lS " (r n+1 r n+2 - ^r n -ir n +i - gr n _ 2 r„_i) . (2) 



Both sides have been multiplied by e~ l8n and we have introduced the new variable S n (t) = Y^j=o® n +j(t)' 
1, . . . , N — 2 with the boundary conditions SN-2(t) — <SV-i(i) = SV(i)- S n is the natural phase parameter because 
the evolution of the model is invariant under any rearrangement of the phases in which the values of S n are conserved 
p2| , ^9|j20| . By separating Eq.|2| into real and imaginary part, one obtains 

f n = -vk 2 n r n + sin(S„) • R n (3) 
6 n = cos(S n ) ■ R n /r n . (4) 



where 



Rn = k n ( r n+1 r n+2 - 7r„_!r„ +1 - ^r„_ 2 r„_i ) • (5) 



R n is the real valued coupling from the nearest shells on the n th shell. Combining Eq.|| and EqJJ one eliminates the 
coupling from the nearest shell and obtains an equation for the time derivative of the phases: 

6 n = C0t(5 n ) + uk n j (6) 

If we ignore the two neighbouring phases in Eq.^ and replace S n with 6 n , simple linear stability analysis of Eq.^ 

gives that the phases are attracted/repelled from the fixed points 9 n — ±-| depending on whether + uk%\ is 

negative/positive. Nearly the same stability conditions is found for S n when the effect of the neighbouring phases 
has been included in Eq.|| Jlfif] . Direct measurements of the stability of the phases show excellent agreement with 

the stability predicted by the sign of + uk^j . The only coupling from the phases on the equation of the moduli 

(Eq.||) is the factor sin(iSVi) which is close to minus one during most of the evolution. If we set sin(S , „) equal to -1, 
the effect of the phases is removed from the equation of the moduli which then becomes a GOY model in terms of 
real variables. By comparing the evolution of the complex and the real valued GOY model we get an estimate of the 
effect of the phases on the evolution of the complex GOY model. It turns out that the phases have roughly no effect 
on the appearance of bursts since the real valued model create approximately the same bursts as the complex model. 
We therefore begin by studying the bursts of the real valued model. The main feature of bursts is the attraction and 
repulsion of the amplitudes to r n — 0, which is a trivial but important fixed point for the model. This dynamics 
appears as a result of a balance between the viscosity term and the coupling term. In the simplest form the real 
valued GOY model can be written as 

f n = V n + C n (7) 

where V n = —vk\r n is the viscosity term and C n = —k n (r n+ ir n+ 2 — \r n -ir n+ i — |r n _2^n-i) is the coupling term, 
and where the forcing is neglected because we focus on the dynamics of the high wave numbers. As V n and r n are 
proportional with opposite signs, Eqj^ can also be written as V n = —vk\{y n + C n ). The values of V n are shown in 
Fig. 3 by a grid of arrows. Only the V„-field is shown by arrows as this determines the change in the attraction (or 
repulsion) to (or from) the fixed point. Furthermore, the values of C n are not universal. Also shown are trajectories 
of (V n ,C n ) during a burst, both of the real valued model (solid) and the complex model (dotted) (for the complex 
model the moduli are drawn). The dashed straight line show where the flow vanishes (f n — 0). First we notice that 
the qualitative similarity between bursts in the real valued and complex model show the weak effect of the phases on 
bursts. From the flow field in Fig. 3 we see that without variations in the coupling term, the amplitude will stabilize 
at the dashed line. 
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Each stage of the dynamics is labelled in Fig. 3. During the attraction, labelled A, the viscosity term and the 
coupling term balance each other with a slight dominance of the viscosity term. Because the trajectory approaches 
the fixed point r n — with decreasing velocity and because \V n \ and \C n \ are much larger than \r n \, the trajectory 
always approach r n = tangentially to, but slightly below, the diagonal V n = —C n . When approaching r n = the 
trajectory can only be kicked away by variations in the coupling term, and the absence of these variations makes the 
amplitude stabilize close to zero. The delicate balance between V n and C n is therefore responsible for the long lasting 
laminar regimes between the bursts. 

As soon as a burst approaches from the lower shells the coupling term becomes large and the trajectory is forced 
away from r n — in a given direction (labelled B in Fig. 3) into a regime of positive f n . This direction is not universal; 
by construction it depends on the value of the amplitudes in the neighboring shells |21|. At the end of the repulsion, 
labelled C, the coupling becomes less dominant over the viscosity, and this is seen in Fig. 3 as a turning towards the 
dashed equilibrium line. As this line is crossed, the amplitude reaches its maximum, after which it is again attracted 
towards zero by the dominance of the viscosity. The motion is highly exited by the burst and the amplitude oscillates 
around the line. The oscillations are damped away and the amplitude settles again close to the dashed line where the 
picture is now repeated again beginning at A. 

From the above scenario it is possible to give an explanation to why bursts are created in the GOY model: The 
amplitude gets trapped at the fixed point in zero and can only be released when a burst arrives from the lower shells. 
As soon as a weak burst is created at the low shells, it continues all the way to the highest shells because the stability 
of higher shells are changed by the approaching burst. The intermittency is created by a "domino" effect through the 
shells. The highest and the lowest shells evolve differently because the effect of the viscosity reduces towards the low 
shells. This reduces the attraction of the amplitudes towards zero, which makes it less possible for bursts to occur 
and gives instead a slow random walk dynamics with a Gaussian statistics. The low shells therefore produce slow 
random perturbations which propagates up through the shells and release bursts at large shells (small scale). 

In conclusion, we have described the mechanism of the creation of intermittent bursts in the GOY model. The 
results show that the creation of a burst is determined by a delicate balance between the viscosity and advection 
terms. We therefore believe that a similar scenario might be present in other intermittent, turbulent systems and also 
in experiments. Our main observation is that a burst is associated with a "fingerprint" : The amplitudes of the high 
wavenumber modes vanish before the burst. An experimental time signal, say from hot wire measurements, might 
indeed show similar characteristica. We are in the process of investigating this using wavelet analysis around the 
bursts. Similar work in this direction has also been done recently by Camussi and Guj ]2^ j. 

We are grateful to P. Bak, T. Bohr, S. Ciliberto, T. Dombre, K. Hansen, J. Kockelkoren, G. Zocchi for discussions 
and suggestions. 
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FIG. 1. The temporal evolution (Eq.jjj) of the logarithmic of the modulus of modes corresponding to the highest shells 
(ris, . . . ,ri9) for a time-span between two bursts. The uppermost curve corresponds to shell 15, while the lowest to shell 19. 
The parameters for the numerical integration are listed in the text. 

FIG. 2. The time average of P„(t) = \0 n — 9 n --s | normalized by n versus the shell number n. Note that the phase organization 
corresponding to P n (t) ~ holds very well for the highest shells. The rise of the graph at n — 19 is due to a boundary effect. 

FIG. 3. Trajectories of V n vs. C n during a burst for the real valued model, Eq. ^ with forcing added (solid curve), and the 
complex model, Eq. jj] (dotted curve). On the same graph is shown the flow field of the viscosity V„ visualized by arrows. The 
dashed line shows where V n = 0. The labels represent the attraction to the fixed point r n — (A), the repulsion away from the 
fixed point (B), and maximum amplitude of the field (C). The arrows on the trajectories indicate the direction of the temporal 
evolution. 
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